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Placing quantum emitters close to a metallic nano-wire, an effective interaction can be achieved over dis- 
tances large compared to the resonance wavelength due to the strong coupling between emitters and the surface 
plasmon modes of the wire. This leads to modified collective decay rates, as well as to Lamb and dipole-dipole 
shifts. In this paper we present a general method for calculating these level shifts which we subsequently apply 
to the case of a pair of atoms coupled to the guided modes of a nano-wire. 
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I. INTRODUCTION 

Coupling quantum systems to a common reservoir results 
in an effective interaction between these systems. In the 
case of spontaneous emission of light from quantum emitters 
this manifests itself in phenomena such as superradiance and 
dipole-dipole interactions. In free space these effects quickly 
disappear as soon as the average distance between the emitters 
exceeds the resonant wavelength. The situation changes dra- 
matically, however, if the most dominant reservoir modes are 
reduced to one or zero spatial dimensions as in the case of a 
nano-wire or a single-mode resonator. Here interactions over 
large distances can emerge. Both the collective decay rates as 
well as the effective interaction Hamiltonian are fully deter- 
mined by the dyadic Green's function of the electromagnetic 
field characterizing the response of the tailored reservoir. Cal- 
culating the collective decay rates requires to determine the 
Green's tensor at a given frequency, typically the resonance 
frequency of the involved emitters. However, in order to deter- 
mine the level shifts one must perform a principal value inte- 
gral over the whole positive frequency spectrum which raises 
serious additional difficulties. 

In this paper we present a method to simplify the calcula- 
tion of the Lamb and dipole-dipole shifts. One of the key steps 
is to introduce imaginary frequencies and extend the integrals 
over frequency into the complex plane. There are also other 
situations when introducing imaginary frequencies and us- 
ing complex frequency integrals helps to simplify the original 
problem, for example in case of calculating Casimir-Polder 
otentials and investigating van der Waals interactions (0], 
, JH] and references therein). We also introduce a special 
Kramers-Kronig relation which, combined with the aforesaid 
integral extension lets us transform the original expression for 
the level shifts into a much more convenient form. 

We apply the method to a particular example where a pair 
of quantum emitters are coupled to the surface plasmon modes 
of a nano-wire. This system is interesting because it enables to 



attain strong atom-field coupling and single-site addressability 
at the same time. Because of the small transverse mode area of 
the plasmons of a metallic cylinder with a sub-wavelength ra- 
dius, a strong Purcell effect arises between the plasmons and 
a single emitter placed close to the wire (0], JH, J^t). The 
effect of strong coupling to the plasmon modes of a waveg- 
uide has been studied for various specific scenarios (01, JU). 
The system of a single emitter coupled to a wire has been pro- 
posed as an efficient single-photon generator (01), as well as 
a single-photon transistor (|9|) and the coupling has been ex- 
perimentally demonstrated ( 1U0II ). Having a pair of emitters 
coupled to the guided modes, we derived an inter-emitter dis- 
tance dependent superradiance effect in Here we calculate 
the dipole-dipole level shifts of the two-atom system. 



II. GENERAL METHOD 

Let us consider a system of N two-level quantum emitters, 
characterized by the lowering and raising operators frj and <rj 
coupled to a common electromagnetic reservoir. Under con- 
ditions that permit the dipole-, rotating wave-, and Markov 
approximation tracing out the reservoir modes leads to a mas- 
ter equation for the atoms 
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The first, Hermitian term of the rhs contains the Suj nn Lamb 
shifts and Scj m7 t n dipole-dipole shifts. In the dissipative term 
we find the single-atom decay rates T nn and the inter-atomic 
decay couplings T m ^ n . Their explicit form is: 
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Here, Gij(f m ,f n ,uj) is the {i,j} component of the Green's 
tensor for the electromagnetic field including the interaction 
with a passive medium such as a nano-wire, evaluated at fre- 
quency u> and at positions r m and r n . It fulfills the Maxwell - 
Helmholtz wave equation 



~e(f,u) G(f t y,w) = I8(?-?), 
cr J 

(4) 

with the proper boundary conditions. e(r, iS) and fi(r, ui) are 
the relative electric permittivity and magnetic permeability, 
respectively. In case of transition frequencies uia for which 
the rotating wave approximation is valid and which are far 
from the ultraviolet domain, we can use the full (vacuum plus 
material part) Green's tensor for calculating T mn . For cal- 
culating 5uj mn one has to perform however an integral over 
the whole spectrum. Since the atom-field coupling is treated 
in a non-relativistic way, this approach does not take into ac- 
count the relativistic high-frequency contributions correctly. 
A well known consequence of this is that the vacuum level 
shifts (Lamb shift) obtained from © are incorrect. Here a 
rather involved calculation based on relativistic quantum field 
theory is required, taking into account all possible transitions 
of the emitter and including proper renormalizations. How- 
ever, if we are interested only in the changes produced by the 
presence of a medium, this problem can be circumvented. The 
medium tailors the reservoir modes only within a certain fre- 
quency range and becomes transparent in the high-frequency 
domain. Thus, calculating the effects relative to the case of 
free-space vacuum, i.e using G med = G — G vac instead of the 
full G, the above equations give correct expressions also for 
the Lamb or dipole-dipole shifts relative to vacuum. In other 
words calculating only the material-induced level shift intro- 
duces an automatic renormalization and lets us get rid of the 
ultraviolet divergences. 

Deriving (f2]i and (0 we also used a Markov approximation. 
This is possible as long as the calculated decay rates T mn (oj) 
and shifts Su) mn (u)) depend only slowly on frequency u>, i.e. 
do not change appreciably over frequency ranges of the order 
of r„ m and Su) mn , It should be kept in mind that even if 
the spectral response of the medium is flat, retardation effects 
can cause a spectral dispersion of the Green's tensor at two 
different positions r m and r n with a characteristic width given 
by c/\r m -r n \ 01]. 

Even though ultraviolet divergencies are eliminated in (f3) 
by considering only the changes due to the medium, the ex- 
pression still contains an integral which is rather difficult to 
calculate. In the following, we will give a general method for 
simplifying this expression. For this we shall use a general- 
ized Kramers-Kronig relation. 



Kramers-Kronig relations (see Appendix[A]) apply, e.g 
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Because the Green's function inherits the symmetry 
G(— (J*) = G* (w) from e(ui), equation © can also be written 
in the form 



Re{G{tu A )} 



Jo - uj 2 a 



(6) 



An important step in the derivation of the Kramers-Kronig re- 
lation is the integration over the semicircle (~r> contribution) 
in the complex upper half plane. As this integration is done 
for large \ui\, using the vacuum Green's tensor ( |[l"2"ll ) 
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is a good approximation. At large \uj\, it goes as 
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To perform the ^ integral we use lu = |w|e M and approach 
the two points where the part joins the real axis (S < a < 
7T — 6), We then find 
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Taking \uj \ to infinity and 8 to zero, we get 
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Because the integrand on any point of the -r^ contour part 
goes exponentially fast to 0, the integral will vanish on this 
part of the contour even if we multiply the integrand with a 
polynomial of u>. In particular we find that the variant of the 
Kramers-Kronig relation 



A. Kramers-Kronig relations for the Green's tensor 

The full Green's tensor G(w) does not have any poles on 
the upper complex half-plane because of causality. Thus 
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holds as well. For the reasons stated above, ((5} and dTTT > - 
being true for G and G vac - are valid for the material contri- 
bution G med also. 
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B. Medium contribution to the Lamb and dipole-dipole shift 

We can rewrite the principal value integral in ((3) as 

j 1=P / du— X *V - (cv + oja), (12) 
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where we replaced the full G by G med = G — G vac , i.e. the 
contribution of the medium. We can now substitute the variant 
of the Kramers-Kronig relation (fTTT i into ( fT2l which yields: 
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Now, we will try to eliminate the principal value from the sec- 
ond term of (fT3l . To do this, we will have to transfer the in- 
tegral from the real axis to the imaginary axis. This, however, 
shall have its advantages because the Green's tensor behaves 
much more smoothly for purely imaginary frequencies: oscil- 
lations become exponentially decreasing functions. 
Resolving the principal value 
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the second term in ( fT3l assumes the form 
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Because the integral goes from to oo on the real axis, we 
can create a closed contour in the upper right quarter of the 
complex plane. The integral over the curved part (C) of the 
contour will again disappear, so we can write 
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Using ui = in on the imaginary axis, we get 
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Substituting ( fT71 ) into ( fT3l and then that into ([3]), we get for 
the Lamb and dipole-dipole shift 
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In this form we no longer have to worry about the principal 
value. As an additional benefit, transferring the integration 
to the imaginary axis makes the Green's tensor much better 
behaved (exponential decay instead of oscillations) which is 
very useful when calculating the shift by numerical means. 



C. Origin of the integral term 

The expression for the dipole-dipole shift found in the lit- 
erature (for example, lfl3l . llTill . IU5ll ) usually involves only 
the first term of the right hand side of (fT8~t or, equivalently, 
the shift is proportional to the real part of the electric field at 
the atomic frequency. This is, however, an approximation that 
relies on the assumption that only frequencies around uja con- 
tribute to the dipole-dipole shift considerably. We can obtain 
this result if we only keep the first two terms in (H4l . arguing 
that the other terms have their greatest contribution at around 
w = —uja which is not contained by the region of integration 
in ( fT2l . However, there is a fundamental reason why in gen- 
eral the integral term in ( fT8l must be present. Starting from 
(fl3l . we can write 
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Resolving the integral term on the rhs of Eq. ([19), we get two 
terms, one of which is identical to the one on the lhs of the 
equation. Thus, one obtains 
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Changing the integration variable on the rhs of Eq. (120) to 
uj = —uj and using the fact that Im{Gy (uj)} is an odd func- 
tion of uj, we get 



50 uj 2 lm{G"' ed (uj)} uj 2 
duj^r 1 lJ 1 n = 7r^-Re{G^ ed (uJA)} 



UJ — UJA 



duj^r- 



lm{G% ed (u>)} 



(21) 



uj — uja 



which, rearranged, gives us the variant of the Kramers-Kronig 
relation (fTTT i 



P / du- 



j 2 lm{G\Y d {uj)} 

: 2 uj — uja 



= Tr^-Rc{G™ ed (uj A )}. (22) 



Thus, one sees that the discussed term in (fT9l is necessary in 
order to fulfill the Kramers-Kronig relations. Note that the 
integral term on the rhs of ( |20l can be transformed into the 
one we got in ( fT8l ), making use of the well-known property of 
G(r*i, ri, w) to be real-valued for purely imaginary frequen- 
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FIG. 1 : Collective atomic decays and dipole-dipole shifts of a pair of 
interacting two-level atoms, described in the Dicke basis. If there is 
a strong effective coupling (high T12) between the atoms, the differ- 
ence between the collective decay rates becomes large, leading to a 
superradiance phenomenon. 



As we will see later on for the set-up discussed here and 
similar situations the second term in ( fT9] l is indeed negligi- 
ble for emitter separations much larger than the characteris- 
tic length (i.e., the wavelength). However, when the emitters 
get close enough, the static (u = 0) contribution in <5wi2 in- 
creases substantially, and we cannot neglect the integral term 
any more. Also there may be situations such as the coupling 
of quantum emitters over macroscopic distances through a 
negative-index material, where this simple rule may not hold. 



III. A PAIR OF EMITTERS NEAR A NANO-WIRE 

In the following, we will apply the method introduced 
above to a particular example. From the interaction of two 
emitters, each coupled to the basic surface-plasmon (SP) 
mode of a metallic wire of sub-wavelength radius, emerges 
a Dicke superradiance effect that depends on the inter-emitter 
distance (la]). As seen in Fig[T] the full collective atomic de- 
cays are Tu ± Ti2 and the singly excited levels get a ±6ui2 
wire-induced dipole-dipole shift. 

Trying to calculate OJ directly introduces difficulties since 
one has to deal with a principal value integral. This is espe- 
cially a problem if we have to perform the integral numerically 
(which is usually the case by non-trivial geometries), because 
we have to know the behaviour of the G med (ui) around u> = 

uj a . As described in @|, G med (uj) = f™dk z G med (u;k z ), 
which is in this case the scattered part of the Green's tensor. 

Although we know the analytic form of G med (k z ;uj), it is a 
rather complicated function and we cannot analytically inte- 
grate it. Using the method described in the previous section, 
however, circumvents these difficulties and lets us perform the 
much simpler integration in ([T8T l where, in addition we have to 
substitute purely imaginary frequencies into the Green's ten- 
sor, by which we get a more well-behaved function. 
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FIG. 2: Dipole-dipole shift S0J12 due to the presence of the wire and 
the Ti2 coupling between the atoms, scaled by the total single atom 
decay rate Tu, as a function of inter-atomic distance A 2 in units 
of the vacuum radiation wavelength Ao. The wire radius is O.OlAo 
and the distance of the atoms from the wire axis is 0.015Ao- &j"f pr 
is a good analytic approximation of <5wi2, for distances comparable 
to Ao and higher. Here, 8LJ12/T11 shows an oscillating behaviour 
shifted by n/2 with respect to T^/Tii. For inter-emitter distances 
comparable to the wire radius, there is a substantial increase in Slj\2 
and it begins to strongly deviate from Sij±% pr . 



We calculated the dipole-dipole shift resulting from the 
presence of the wire using the full scattered Green's tensor. 
We compared it to an analytic approximation used by ifTill . 
based on a single-plasmon resonance model, a derivation of 
which is given in AppendixlBl 
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^(w A ) 7 (w A )e^ A *sin(*J , Az). (23) 



Here A and 7 are the amplitude and width of the Lorentzian 
fit to the plasmonic resonances (see Appendix [BJ and k v z l is 
the longitudinal component of the wave vector of the plasmon 
mode. Fig. [2] shows the results of the calculations. For inter- 
emitter distances larger than the vacuum radiaton wavelength, 
SuJi2 Pr i s a good approximation to the exact wire-induced 
dipole-dipole shift. It only deviates from the sinusoidal be- 
haviour when the inter-emitter distance becomes comparable 
to the wire radius: in this case, the atoms begin to feel the 
three dimensional nature of the wire and the quasi- ID cou- 
pling approximation (i.e., (5w°f pr ) is no longer valid. How- 
ever, because the wire is quite thin, this typically happens at 
distances well below the vacuum radiation wavelength which 
means that in this regime the emitters are already strongly 
interacting through the vacuum as well. So we can safely 
say that the analytic approximation works well if the inter- 
emitter distance is above the vacuum radiation wavelength of 
the emitters. In the regime where the exact calculations are 
well approximated by 5w°2 Pr , Sui^/Tn oscillates with the 
same period as r 12 /r 11 only with an additional tt/2 relative 
phase shift fl4Tl . This means that for inter-emitter distances 
yielding extrema for T^/Tn, i.e., where the symmetric or 
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FIG. 3: The first term and the second, integral term of the dipole- 
dipole shift, as seen in i I St and the actual wire-induced dipole-dipole 
shift Suj\2, scaled by the total single-emitter decay rate Tn . The wire 
radius is O.OlAo and the distance of the atoms from the wire axis is 
0.015A(). Because of the increasing static (ui = 0) contributions in 
the Green's tensor, for small inter-emitter distances the integral term 
becomes comparable to the first term, becoming no longer negligible. 



antisymmetric transition is superradiant, IS") and \AS) are de- 
generate. On the other hand, when T^/Tn = 0, this degen- 
eracy is lifted by \8lu\2 1 being maximal. The extrema of \8lu\2 1 
are 0.5Fn at most. The decay of the amplitude of the oscil- 
lations for both 8UJ12/T11 and r^/Tii is caused by ohmic 
losses in the metal, represented by 7. Thus, the interaction 
always makes a distinction between the symmetric and the 
antisymmetric transition: either by the different decay rates, 
or by the lifted degeneracy of 15} and \AS). 

As seen in Fig[3j the closer the emitters are to each other, 
the more substantial the integral term of ([18} becomes. This is 
in accordance with the arguments made in Sec. Ill CI namely 
that the decrease of the inter-emitter distance enhances the 
static contribution of the Green's tensor. For small enough 
distances, the integral term is not negligible any more. 



lated. We discuss the appearance and importance of an inte- 
gral term in the derived expression, that sweeps across purely 
imaginary frequencies, and is usually neglected in the litera- 
ture. We apply the method for calculating the dipole-dipole 
shift of a pair of atoms coupled to the guided surface plas- 
mon modes of a metallic nano-wire and compare it to a quasi- 
1D analytic approximation. The results show that for inter- 
emitter distances comparable to the wire radius 5lu\2 becomes 
considerably larger than the single-atom decay rate (Tn) and 
the approximation doesn't hold. However, for larger distances 
the shift is very well approximated by the analytic expression, 
and <5wi2/Fii shows an oscillatory behaviour having roughly 
half the amplitude and the same period as r^/rn, as well as 
an additional phase shift of ir/2. Thus, the interaction always 
makes a distinction between the symmetric and antisymmetric 
transition of the 2-atom system, either by the modified collec- 
tive decay rates (superradiance) or the lifted degeneracy of \S) 
and I AS*}. We also take a look at the behaviour of the integral 
term mentioned above and conclude that it indeed becomes 
substantial for distances comparable to the wire radius. 
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Appendix 

Appendix A: Kramers-Kronig relation generally 

First, let us look at the Kramers-Kronig relations in case of 
a general, complex-valued function f(oj) which is analytic in 
the upper complex half-plane. According to Cauchy's theo- 
rem, 



f(u A ) = limi 

(5->0+ 27T1 



(Lu- 



lua — i<5 



(Al) 



IV. SUMMARY 

In the present paper we discussed the effects of a tailored 
reservoir on the Lamb and dipole-dipole shifts of quantum 
emitters coupled to a common radiative reservoir. By con- 
sidering only shifts relative to those in vacuum all complica- 
tions arizing from ultra-violet divergencies and off-resonant 
contributions from other transitions were eliminated, reduc- 
ing the problem to the calculation of an integral of the elec- 
tromagnetic Green tensor over all frequencies. We presented a 
method that greatly simplifies this calculation by transforming 
the original expression containing a principal value integral 
over frequency into an ordinary integration over the imaginary 
axis. The method does not imply any specific configuration or 
system, so it can be used in a wide variety of problems where 
level shifts due to dipole-dipole interaction have to be calcu- 



and the contour (containing the point lua + iS within) is closed 
in the upper complex half-plane. If we assume that the path 
integral of /(lu) is nonzero only over the real axis and disap- 
pears on the other parts of the contour (that is, on the complex 
plane) if we extend it to infinity, we can write the integral as 
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And so: 



i r°° 

f(u A ) = -P / 

7T1 J-oo 

or with real and imaginary parts 
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(A6) 



because in this case the only substantial channel that couples 
the emitters are the surface plasmon modes. 
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Because of the translational invariance of the wire in the z 
direction, G m is symmetric in k z . Thus, substituting in (IB2b 
we get 



Appendix B: Analytic approximation of the dipole-dipole shift 
for two emitters coupled by a nano-wire 



Im{G™ ed (r u r 2 ,uj)} = Im{2iriAje ik *' Az e-~< Az }. (B4) 



The Green's tensor for an infinitely long, cylindrical wire 
can be calculated as given in |@]. For the scattered part we can 
formally write 



dk z G med {r u r 2 ,u;k z ) (Bl) 

oo 



where we know Q med analytically. The atom-wire coupling 
is strongest if the dipole moment of the emitters point in the 
radial direction. If the cylindrical coordinates of two emitters 
only differ in their z component (where z is the symmetry axis 
of the wire) one finds 



G™ d {fi,f 2 ,uj) 



dk z e lk ^ z GZ ed {run^-k z ) (B2) 



where Az = \z 2 — Zi\. In case of a thin wire (radius well 
below the vacuum radiation wavelength of the emitter), and 
small emitter-wire distance (in the order of magnitude of 
the radius), the plasmonic contribution becomes dominant in 
G™ ed (fi,fi ) (J] k z ). In this case, the imaginary part of the 
Green's tensor is well approximated by two Lorentzian fits, 
centered at ±fc^', i.e., the z component of the wave vector of 
the plasmonic mode. This approximation is valid for inter- 
emitter distances larger than the vacuum radiation wavelength 



With this, we can now express the wire-induced single-emitter 
decay rate and emitter-emitter coupling, respectively 
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l -A{Lu A )l(uj A )e-^ Az cos(fcP z Az),(B6) 

as well as the wire-induced dipole-dipole shift, according to 
(HI: 
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According to the discussion in the paper, for inter-emitter dis- 
tances larger than the vacuum radiation wavelength the inte- 
gral term in (|B7t can be neglected, so in the end we arrive to 
the analytic approximation 

§ ^ pr ^ _^ i(Wj()7(wA)e - 7 A Zsi4 ,l A4 (Bg) 
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